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The three-dimensional integer- valued lattice gauge theory, which is also known as a "frozen super- 
conductor," can be obtained as a certain limit of the Ginzburg-Landau theory of superconductivity, 
and is believed to be in the same universality class. It is also exactly dual to the three-dimensional 
XY model. We use this duality to demonstrate the practicality of recently developed methods for 
studying topological defects, and investigate the critical behaviour of the phase transition using 
numerical Monte Carlo simulations of both theories. On the gauge theory side, we concentrate on 
the vortex tension and the penetration depth, which map onto the correlation lengths of the order 
parameter and the Noether current in the XY model, respectively. We show how these quantities 
behave near the critical point, and that the penetration depth exhibits critical scaling only very 
close to the transition point. This may explain the failure of superconductor experiments to see the 
inverted XY model scaling. 



PACS numbers: 74.60.-w, 64.60.Cn, 64.60.Fr, 11. 15. Ha 

I. INTRODUCTION 

Duality arguments suggest that the superconductor- 
insulator phase transition should be in the same univer- 
sality class as the three-dimensional XY model, but with 
an inverted temperature axisi2ii^ This means that the 
superconducting phase maps onto the symmetric phase, 
and the normal phase onto the broken phase. Further- 
more, the dual counterpart of the order parameter of the 
XY model is a non-perturbative field that creates the 
Abrikosov-Nielsen-Olesen vortices in the superconductor. 

The most direct prediction of the duality is that the 
critical magnetic field H c \ , or equivalently the vortex ten- 
sion T, should scale with the XY model critical exponent 
^xy ~ 0.6723, - In practice, it is easier to measure the 
scaling exponent v' of the penetration depth A, but for 
that quantity, the theoretical picture is less clear. Both 
v' « 0.33£ and v' « O.S 4 *^ were suggested before the 
prediction eventually converged to v 1 = ^xy- 8 

Ironically, two different experiments with 
YBa2Cu307_,5 high-temperature superconductor have 
produced results that are each compatible with one 
of the earlier suggestions, namely v' — 0.34(1)^ and 
v' = 0.45 . . . 0.5ia 

Similarly, it has turned out to be difficult to confirm 
the duality in Monte Carlo simulations of the Ginzburg- 
Landau theory away from the London limit^ as they 
also seem to favour v' sa 0.3. In the London limit, where 
the duality is on a firmer footingji2*i£ simulationsi^i^ 
give v 1 w 0.67 through indirect measurements, though. 

There are two principal reasons for all this confusion: 

1 . The duality relation is expected to apply quantita- 
tively only in a narrow temperature interval near 
the critical point, where the vortex tension/ fc^T 



is well below the inverse correlation lengths of the 
scalar and photon fields. 

2. The duality relates the fundamental fields and 
thermodynamical densities of one theory to non- 
perturbative and non-local objects in the other. 
These are difficult to measure both in experiments 
and in numerical simulations. 

In the London limit, the difficulty (1) is alleviated, 
as the scalar correlation length is very short, but in 
this paper we go even further and take another limit 
to obtain the Abelian integer-valued lattice gauge the- 
ory, "frozen superconductor" (FZS) X For this theory, the 
duality transformation can be carried out exactly, and 
yields precisely the three-dimensional XY model with the 
Villain actionii2ii§i±£ The duality is therefore exactly 
valid at all temperatures. 

The FZS does not have a scalar (Higgs) field to drive 
the transition, but it still has a transition between a low- 
temperature superconducting phase with massive gauge 
field excitations, and a high-temperature massless phase. 
(Microscopically, this transition is due to the "freezing" 
of the discrete gauge variables below the critical temper- 
ature; hence the name frozen superconductor.) This is 
mapped to the symmetry breaking transition of the XY- 
Villain model, with inverted temperature (T <-> 1/T). 
Since the duality is valid at all temperatures, the prob- 
lem (1) above is completely avoided. 

In this paper we study the duality relation in detail 
with lattice Monte Carlo simulations. Our aim is not 
to numerically verify the duality — it is, after all, a 
mathematical identity. Rather, the purpose is to iden- 
tify "good" observables and possible pitfalls in Monte 
Carlo studies of the duality relations of this type, espe- 
cially bearing in mind the problems encountered in the 
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GL theory simulations. By this we mean, first of all, that 
we want to construct pairs of observables in the frozen 
superconductor and in the XY model which are exactly 
dual to each other and which are sensitive to the critical 
behaviour. The observable in the frozen superconductor 
must be defined also in the full GL model. It is, however, 
equally important that the observables can be measured 
to a high accuracy in both theories. The exact duality of 
the models we consider allows us to check this and iden- 
tify the errors caused by finite-size effects and possible 
inefficiencies caused by the numerical methods used. In 
this way, we hope to find out how to avoid the difficulty 
(2) in more realistic cases. 

The main body of the paper is concerned with the cal- 
culation of gauge invariant (non-local) order parameters 
of the FZS; the vortex tension T, the photon mass m 7 
and the magnetic permeability Xm- We relate these ob- 
servables to their duals in the XY model, and study the 
critical behaviour. 

This paper is organized as follows. Section 2 contains 
the definition of the two models and discusses the duality 
between them. Section 3 introduces the observables and 
Section 4 discusses their qualitative behaviour in the two 
phases. In Section 5, the location of critical point is de- 
termined using the spin correlator in the XY model, and 
Section 6 deals with the definition and the measurement 
of its dual counterpart, the vortex tension. In Section 8 
we study the susceptibilities, i.e., the photon correlator 
in the momentum space, and discuss the mass determina- 
tion from them, and in Section 7, we extract the photon 
mass directly from the exponential decay of the correla- 
tion function in the coordinate space. Finally, Section 9 
lists the main findings and concludes the paper. 



II. DUALITY 

We start by formulating the three-dimensional 
Ginzburg-Landau theory in the London limit on a three- 
dimensional cubic lattice. The theory consists of a real 
valued gauge field A x ^ defined on links and spin angles 
9 X defined on lattice sites. The partition function is 



where 



£gl,x = - \ E F la + K E s - *- ^) > ( 2 ) 
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and 
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(3) 



and s(x) is a periodic function with period 2tt and 
minimum at i = 0. The standard choice for s(x) is 
s(x) = — cos(a;), but here we shall use the Villain formAS. 



s(x) 
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-In y \ exp y — -{x — 2irk)' 



(4) 



In this paper we shall consider two limits of this theory. 
The frozen superconductor (FZS)^i^ which we also refer 
to as simply "the gauge theory," is obtained by taking 
k — > oo and defining (3 — 4ir 2 /q 2 . This leads to the 
partition function 



z F zsW = E ex p I -f E R - 
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X,'1>J 



where 



P — T 



(5) 



(6) 



and the link variables = qA Xi i/2TT take integer values. 
This model has two phases: the Coulomb phase ((3 < (3 C ) 
and the superconducting phase (/3 > (3 C ). 

We also consider the limit q — > 0, in which we recover 
the three-dimensional XY model. The gauge field A Xj i 
decouples, and the non-trivial spin part of the partition 
function becomes 



Zxy(k) = J £>0exp ^-KE s ( 0x +*-^)j 



(7) 
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This model has a broken phase at k > k c and a symmetric 
(1) phase at k < k c . 



We shall now show that the partition functions © and J7J) are dual to each other, i.e., proportional to each other 
if P = 1/k. Introducing a real vector field /i x ,ir^ we can write Eq. J7J) as 



Z xy (k)oc / mo/i, E ex P 

J l. 



E 



— h 2 - ih A XY 

ZK 



(8) 



where we have introduced the Noether current of the XY model 



(9) 
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The integration over 9 yields a delta function S(V ■ h), where we have defined the lattice divergence 

V • h x = 2J (hx,i - hx-i,i) ■ 



(10) 



The summation over fc Xj i restricts /i x ,i to integer values, and we obtai: 
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Zxy(k) oc *v-h,o exp ( ^ 



(ii) 



which is the partition function of an integer-valued and sourceless vector field. In an infinite volume, we can interpret 
the vector field h x i as the integer valued flux through the dual lattice plaquette pierced by link (x, i), and write 



(12) 



with -Pxjfc as in Eq. 0. 

Thus, identifying (3 — 1/k, we recover the partition function in Eq. 1,12,17 This shows that the two limits of 
Eq. JTJ), k — > oo and q — > 0, are dual to each other. In particular, the duality relation implies that the gauge theory 
has a phase transition of the XY model universality class at f3 c — 1/k c , and the Coulomb and superconducting phases 
are mapped to the broken and symmetric phases, respectively. Note, however, that duality only maps the two limits 
onto each other; the Ginzburg-Landau theory is not self-dual for finite parameter values. 

The nature of the duality transformation becomes more transparent when we introduce an external, real-valued 



vector field a x ,i in the XY model 



Zxy {k; {a x ,i}) = J DO ^ exp 



Introducing h x ,i as in Eq. (jSJ), but defining h x i = h Xi i — a Xj i, we obtain 



(13) 



Z X y{k; {a X:i }) cx / DODhj exp 



i/t x .,,A*y 



(14) 



Now, as before, we integrate over (9 X , which constrains V • h x = 0, sum over fex,i> which makes /i Xi i an integer, and 
express ft, Xi j in terms of P x .ij as in Eq. 1)12(1. This gives us 



Zxy («; = z fzs (1/k; {a Xii }) , 



where 



Zfzs (P; {a x ,i}) = 2J ex P 



2 ^ V 2 



(15) 



(16) 



Note that in this expression the vector field a is defined on the dual lattice, i.e. a Xj i lives on the link that pierces the 
plaquette P x ,jk- Eq. (|15fl is the basic duality equation, which can be used to relate the observables of the two models 
to each other. 



III. OBSERVABLES 
A. Spin-spin correlator 



Using Eq. (|13|) , we can write this as 

Zxy (k; {a Xti }) 



6' 



Zxy{k) 



(18) 



The basic observable in the XY model is the spin-spin where a X;! : is an otherwise arbitrary fixed integer-valued 
correlation function vector field, but it satisfies the condition 



G x 



(exp[i(0 Xl -0 X2 )]> 



XY ' 



(17) 



V • a = <5 x , Xl - 5 x ,x 2 , 



(19) 
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i.e. it has a source and a "sink" at points Xi and X2, 
respectively. 

The duality relation in Eq. (|T5|) implies that Eq. l (T%|l 
must be equal to 



6' 



eXp I ~ 9 E J- + a 




FZS 

(20) 

The simplest choice for a x ,-t vanishes everywhere except 
on the shortest path of links that leads from xi to X2 , on 
which it has the value of unity. As we shall see later in 
Sec. IVII Eq. H2Ufl then has a natural interpretation as a 
vortex correlation function in the gauge theoryiifiaSi In 
Table [l]we summarize the duality between basic observ- 
ables in the XY model and in the frozen superconductor. 

In the symmetric phase of the XY model, i.e., when 
k < k c , the spin-spin correlator decays exponentially as 



cx e 



<2| 



(21) 



where we call the decay rate m the scalar mass, in accor- 
dance with a field theory picture. By definition, the cor- 
relation length is given by its inverse £ = l/m. Eq. i|20|) 
implies that under the duality transformation, the scalar 
mass becomes the vortex tension defined in Sec. IVII 
T = m. 

In the broken phase of the XY model, where k > K c , 
the spin-spin correlator l|17|l approaches a constant 



lim G XliX2 = M 

|xi-x 2 |— >oo 



(22) 



where M is the magnetization. 

More generally one can also consider higher n-point 
functions in the XY model, which correspond to external 
fields a x j with more sources and sinks. 

We would like to emphasize the difference between our 
approach and the attempts to describe the phase transi- 
tion of the XY model as vortex percolationi2ii22iS The 
"line tension" discussed in that context is defined us- 
ing the length distribution of vortices in the XY model 
and suffers from certain ambiguities, and there is numer- 
ical evidence that the percolation point at which the line 
tension vanishes does not even coincide with the ther- 
modynamic critical point re c ^» In contrast, our T is the 
vortex tension in the gauge theory, is a well-defined, un- 
ambiguous observable and, due to the exact nature of 
the duality, reflects the true thermodynamic properties 
of the system. 



9{x, y, z) + 89. It is convenient to define a periodic vari- 
able 

§(x, y, z) = 9(x, y, z) - — 80 = 0(x, y, z) - zj 3 , (23) 

iV 3 

where js = 89 /N is the average Noether current created 
by the twist. We can then write the twisted partition 
function as 



Zxy = J D9 ^2 ex P 



{**:,*} 



E (? 



+ 3i 



where j = (0,0, j 3 )- 

The helicity modulus is defined as 



T 



N 3 ( d 2 F \ 



1 (d 2 F 



N 1 N 2 {dS9 2 J se=0 V\dfJ h=0 



(24) 



(25) 



where F = — In Zxy is the free energy of the system. It is 
an order parameter, with finite, non-zero values T > in 
the broken phase (k > k c ) and a vanishing value T = 
in the symmetric phase (n < k c ). If the XY model is 
viewed as a model for a superfluid, T is proportional to 
the superfluid density, T cx p s . We can also write T as 



1 1 d 2 Z XY 



VZ XY 
k(1- kx) 



K 

V 



V 



(26) 

where \ = J3 x (Aq ^A^g) is a susceptibility related to 
the U(l) current density of the XY model. 

It is interesting to ask what the helicity modulus T 
corresponds to in the gauge theory. Formally, it is related 
to Eq. (|13|l with a constant imaginary field on — —inji, 



Z 



XY 



(K,jt) = Zxy («; {on}) exp (~Vj? ) , (27) 



and therefore 

- 2 i d 2 z 



T = 



K 



XY 



1 



1 



^FZS 



V Zxy dai. (5 2 V Z FZS 

Using Eq. (JSJ, this is nothing but 




da 2 



Xr, 



(29) 



where we have defined the magnetic permeability Xm as 
the susceptibility associated with the magnetic flux in 
the FZS. 



B. Helicity modulus 

Another important quantity in the XY model is the 
helicity modulus^ which characterizes the response of 
the free energy to a twist of the spins by an amount 89 
along, say, the z direction. Instead of period boundary 
conditions for 9, we would instead have 9(x, y, z + N) — 



C. Photon mass 

From the FZS point of view, the most natural observ- 
able is the photon correlator, defined as the correlation 
function of plaquettes, 

r( Xl .i)( X2j ) = (P xl ,kiP X2 , mn) fzs- (ikl and jmn cyclic) 

(30) 
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XY model 


FZS 


spin-spin correlator 




Eq. (HOJ 


vortex correlator 


scalar mass 


m = l/C 


T 


vortex tension 


helicity modulus 


T 




magnetic permeability 


current-current correlator 




r(x 1 ,i)(x 2 ,j) 


photon correlator 



TABLE I: Comparison of the basic observables in the XY model and in the frozen superconductor (FZS). 



In the superconducting phase ((3 > /3 C ), this correlator 
decays exponentially, and we call the decay rate the pho- 
ton mass m-y. The penetration depth A is defined as 
the inverse of the photon mass, A = 1/to 7 . In the non- 
superconducting Coulomb phase (f3 < j3 c ), F( Xl i)( X2 j\ 
has a power-law decay, which corresponds to vanish- 
ing photon mass m 7 = 0, or infinite penetration depth 
A = oo. 

We write the photon correlation function as 



1 



(xi,i)(x 2 j) 



d 2 Z FZS (ft {a Xii }) 



f3 2 Z F z S ((3) da Xl:i da X2ij 



(31) 



a=0 



which is simply a generalization of Eq. (|28|l to a non- 
constant a x ,i- Using the duality l|15|) . we find 



r 



1 



1 



(x,i)(y,i) - -gS^ySij - (A x jA y J) XY . (32) 

In practice, it is more convenient to study correlations 
between planes rather than individual plaquettes. To 
this end, we interpret the z direction as "time" and label 
it with r. We consider the photon correlation function 
with a non-zero "spatial" momentum p, where we use 
the underline to indicate a two- vector in the (1,2) plane. 
(The zero-momentum correlation function vanishes iden- 
tically.) The correlation function is defined as 



r(r,p) = ^Re J ]T^(P x ,i 2 P y ,i 2 ) 
3 [t x,y 



-ip-(x-y) 



where pi — 2itki/Ni, (i = 1,2) and ki are integers, x = 
(r ,x) and y = (r +r,y). 

It is also useful to consider the three-dimensional 
Fourier transform of the photon correlator 



r y -(p) = £ e ~ ip ' xr V 



(34) 



Eq. I|32|) then translates into 



r«(p) = ~ki - j 2 <Af Y (-p)Af Y (p)) XY . (35) 



Note that 



Xm = lim r 33 (p) 



(36) 



More generally, we can see that there is a direct corre- 
spondence between the photon in the gauge theory and 
the Noether current density A X Y in the XY model. 



IV. PHASE STRUCTURE 

A. Symmetric/superconducting phase 

In the XY model, the symmetric phase (k < k c ) is 
characterized by vanishing magnetization M , and a finite 
correlation length £ = 1/m. The helicity modulus T is 
zero, and therefore Eq. (|26|l implies that x — V K - When 
the critical point is approached, the correlation length 
£ diverges as £ ~ \n — n c \~ u ■ Numerical studies [using 
the cosine action rather than Eq. (J2J] have shown that 
v w 0.6723£ 

In the FZS, this corresponds to the superconducting 
phase (f3 > (3 C ). The above implies that the vortex ten- 
sion T, which is equal to the scalar mass m = l/£ of the 
XY model, is non-zero, and vanishes as T ~ \0 — (3 C \ V 
at the transition point. The magnetic permeability Xm, 
which is equal to T vanishes. However, we can define the 
gauge field susceptibility \A by [cf. Eq. 



XA 



lim 



r 33 (p) 



(37) 



where the underline indicates that p is a two- vector in the 
(1,2) plane. This quantity diverges as the critical point 
is approached, and we parameterize this divergence by 
the exponent ja, 



Xa~\[3-Pc\- 



(38) 



It was argued in Ref. [Tj that 7,4 = v. 

The photon correlator T^ Xl i ^ X2 j^ decays exponentially 
in this phase, and the penetration depth A is therefore fi- 
nite. It diverges at the transition point as A ~ \(3 — (3 C \ V . 
By duality, the current-current correlation length of the 
XY model has the same behaviour. There has been a lot 
of debate in the literature about the value of v' . Orig- 
inally, it was believed that v' — but Kiometzis et 
ali^ later argued that the penetration depth does not get 
renormalized and would therefore have the mean-field ex- 
ponent v' = 1/2. Finally, Herbut and Tesanovio& found 
v' = v using renormalization group arguments, and this 
value was later confirmed in Refs. Il4lll5ll26l with differ- 
ent approaches. However, all these results rely on some 
analytical approximations. They also disagree with the 
results of superconductor experimentsAiS 
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XY: symmetric K < k c 
FZS: superconducting /3 > j3 c 


critical 


K, — K,c 


XY: broken n > k c 
FZS: Coulomb (3 < /3 C 


£, = T' 1 ~ |/t — k c \~" ^ 


G x ,y ~ 


| x _ y |-a+.) 


T = 


A = m^ 1 ~ \k — k c \~" 


r, 3 (p) r 


- Ipr 


A = m^ 1 = co 


X A ~ \P-P C \- 1A 






T ~ Xm ~ |« — Kc 1" 



TABLE II: The behaviour of certain observables in the two phases of the models. 



B. Broken/ Coulomb phase 

From the XY model point of view, the characteristic 
property of the broken phase (n > k c ) is non-zero magne- 
tization M =/= 0. Another signal for symmetry breakdown 
is a non-zero value of the helicity modulus T. Through 
the duality, this corresponds to non-zero magnetic per- 
meability Xm in the Coulomb phase of the FZS. 

In the extreme high-temperature limit j3 — > 0, the FZS 
approaches free non-compact electrodynamics, which is 
given by Eq. 0} in the limit k — * 0. In this case, it is 
easy to calculate \m from the path integral, and we find 

Xm~V/3 when/3^0. (39) 

Using Eqs. and J2§J, we find T — > k and x — ► 0, as 

k — > oo. 

As the critical point is approached Xm vanishes, and 
following Ref. l25l we parameterize this with the critical 
exponent v, 

Xm~\P-(3c\ V , PSPc- (40) 

By duality, the helicity modulus T must behave in the 
same way. It has been argued in Refs. Ill2l27l that v = 
v. The critical exponents are summarized in Table ITT1 

C. Critical point 

At the transition point, the spin-spin correlator in 
the XY model has a power-law decay G XliX2 ~ |xi — 
X2I ~( 1+v \ where the anomalous dimension rj has been 
measured to be 77 « 0.038 & 

Similarly, one can define the anomalous dimension t\a 
for the photon correlator in the FZS by 

ry(p)~|pr, when|p|^0. (41) 

Earlier studies have shown that t\a ~ \^^^^ 

V. LOCATING THE CRITICAL POINT 



K(N) 



.3334 




333 1 1 1 1 1 1 1 1 1 

20 40 60 80 

N 

FIG. 1: Pseudocritical hopping parameter values k(N) as 
a function of the system size N in the XY model Eq. J7J, 
together with a fit to the scaling ansatz in Eq. 14311 . 
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FIG. 2: The scalar mass m in the XY model as a function of 
the dual coupling (3 = 1/k. The solid curve shows the power- 
law fit in Eq. 1441 , and the dashed line the vortex tension fit 
in the FZS [see Eq. <lrJ21 1. 



In order to explore the manifestations of the duality, 
we study both the FZS Eq. © and the XY model Eq. Q 
using Monte Carlo simulations. For the FZS, the simula- 
tion algorithm consists of a single hit Metropolis update 
of the integer link variables, and for the XY model we 
use a Swendsen-Wang type reflection cluster algorithm. 29 



One update sweep for the FZS consists of N1N2N3 single 
link Metropolis hits, and for the XY model out of one full 
Swendsen-Wang cluster update. The number of Monte 
Carlo update sweeps typically ranges from 50000 to 10 7 , 
depending on lattice sizes and coupling constant values. 
The errors are determined by jackknife error calculation. 
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Because of the duality, the critical points of the theo- 
ries are related by (3 C = 1/k c . We determined the critical 
coupling values 



= 0.333068(7), (3 C = 3.00239(6) 



(42) 



by a numerical simulation of the XY model on cubic 
N = Ni = N 2 = N 3 lattices with N = 8-96. For this, 
we used the Binder cumulant method^ Matching the 
fourth-order cumulants of the magnetization measured 
from lattices of size N and 2N gives us an estimate of the 
iV-dependent pseudocritical hopping parameter k c (N). 
The results are shown in Fig. ^ Using the finite size 
scaling ansatz 



K C (N) = K c + 



A K 



(43) 



with u> w 0.8 denoting the subleading XY model scaling 
exponent** we arrive at the numerical result in Eq. l|42Jl . 
The uncertainty in this value is small enough to have 
a negligible effect for the determination of the critical 
indices. 

Most numerical studies of spin models with XY 
universality use the cosine action s(x) = — cos(ai) rather 
than Eq. and subsequently there is only limited ex- 
perience with the Villain action. Nevertheless, there is 
no doubt that the phase transition with the Villain ac- 
tion is of second order and in the same universality class 
as with the cosine action. As an illustration we show the 
scalar mass m — 1/f in the symmetric phase (k < k c ), as 
a function of the "temperature" f3 — in Fig. The 
solid curve in the figure shows a x 2 -fit with the scaling 
law 



m(J3) = - f3 c y 



(44) 



and with a fit with x 2 /d.o.f. = 0.35 we obtain the param- 
eter values — 1.32(4) and v = 0.66(2), in agreement 
with the cosine action results^ The amplitude value 
is important for the comparison with the vortex tension 
in the gauge theory, to be discussed below. 



VI. VORTEX TENSION 

In order to illustrate the numerical consequences of 
the duality we measured the nearest-neighbour spin-spin 
correlation function, 



6' 



= cos 



0) 



XY 



(45) 



in the XY model at (3 = 1/k = 3.02. In the gauge theory, 
this is mapped to the expectation value of 



FZS- 



(46) 



In Fig. we compare both measurements as functions of 
the system size N on cubic boxes with periodic boundary 
conditions. It is evident that Eq. I|12|) . and thereby the 



<cos(e .-e .)> 
— ^ — ±— 




15 20 25 30 N 



FIG. 3: A numerical check of duality relations for the nearest- 
neighbour correlator in Eq. 1451 on a finite lattice as a function 
of the linear lattice size N. 



duality, only is valid for infinite systems, i.e., the bound- 
ary conditions of finite systems do not respect our duality 
arguments. The curves in the figure assume finite volume 
corrections of the form exp(— cN) and extrapolate to the 
same value in the thermodynamic limit. 

Let us now consider the long-distance behaviour of the 
correlation function G XljX2 defined in Eq. 1|17|) . The du- 
ality maps this correlator into Eq. (|20|l in the FZS. In 
principle, we could measure Eq. (|20fl directly, but it turns 
out to be more convenient to use an indirect approach, 
which was introduced in Ref. |2(1 

As we shall now show, Eq. H20[) corresponds exactly to 
a vortex, or an Abrikosov flux tube, in the FZS. Thus, the 
spin-spin correlation length £ in the XY model is exactly 
the tension 7~, the vortex free energy per unit length, 
in the FZS. These flux tubes exist because a supercon- 
ductor resists applied external magnetic fields. For field 
strengths above the critical field H c \ magnetic flux pen- 
etrates the material, and the flux arranges into vortex 
lines, Abrikosov flux tubes. 

In practice, we define the vortex tension as follows: Us- 
ing appropriate modifications of periodic boundary con- 
ditions, to be described below, we constrain the net num- 
ber of vortex lines ny winding around the finite volume 
to, say, z direction. The vortex tension T is then defined 
by 



r 



J_, r Z FZ s{ny = 1) 
N 3 n[ Z FZS {n v = Qy 



(47) 



in the ./V3 — > 00 limit. 

Let us now discuss how the appropriate boundary con- 
ditions are constructed. Labelling the cartesian coordi- 
nates x = (fii,n2,«3) i n-i — l,...,Ni i = 1,...,3, and 
using periodic boundary conditions to all directions 

I(n 1 =Ni + l, n 2 , n 3 ),i = I(n 1 = l.n 2 ,n 3 ),ii 
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W(z) 




FIG. 4: The function W(z) [see Eqs. and EH] measured 
on a 16 3 lattice, with /3 varying from 3.0 (bottom curve) to 
4.0 (top). 



^(ni,ri2,n3=iV3+l),i -^(rii,ri2,r!3 = l),ii 



(48) 



the total magnetic flux $ = 27rny through any planar 
2-dimensional cross section fL » of the lattice 



(49) 



equals zero because of the periodicity. However, if we 



choose a set of links I, 



(ni ,B2 ,ns) 



2 , with ni = 1, n,2 = 1 and 



ri3 = 1 . . . A/3, and use a "twisted" boundary condition 

^(JVi+i,i,n s ),2 = -f(i,i,7i 3 ),2 + nv, n 3 = 1, N 3 , (50) 

the flux through N3 cross sectional areas Ox 2 acquires 
the value $(^1.2) = 2irnv- At ny — 1 this corresponds 
to a single vortex of length L = N3. The vortex line 
forms a closed loop through the z-direction of the box. 

It should be noted that the modification of the bound- 
ary condition does not give the modified link a special 
status: by a suitable redefinition of the link variables Ij, 
the "twist" can be moved to arbitrary location along the 
(l,2)-plane, also away from the boundary. 

In practice, a more convenient way to implement the 
fixed flux is to consider a system with fully periodic 
boundary conditions but with a modified action: Let us 



define a plaquette field which is non-zero only on one 
fixed "stack" of (l,2)-plaquettes, 



if xi = X2 = 1; i = 3, 
otherwise 



(51) 



Now we can define the partition function 



^FZS 



(P,nv) = 



exp 



x . i 



nymy. 



G52) 



which is equivalent to the partition function with the 
unmodified action but with the boundary condition Eq. 
l|5Up. (Again, the choice X\ = x% = 1 here is arbitrary.) 

Comparison of Eq. I|52|l with Eq. (|20|l shows that the 
vortex tension in Eq. I|47() is closely related to the XY 
model spin correlator G Xl . X2 , provided that we identify 
m x ,j with a x ,i- Indeed, the only difference is that m Xj j 
does not have start or end points, but stretches across 
the whole system. However, on a periodic lattice, we 
can imagine moving the source xi of a X) , through the 
boundary to x 2 so that it cancels the sink. This leads to 
a sourceless field which has one field line passing through 
the lattice, and this shows that in the infinite volume 
limit the XY model correlation length £ and the vortex 
tension T are related by 



T =\ji = to. 



(53) 



As an aside, we note that this procedure does not 
make sense for the XY model spin correlators, since now 
G x 

x+e 3 iV3 — Gx.x — 1, again showing that the periodic 
boundary conditions in a finite volume do not respect the 
duality relation. 

In practice, all measurements arc affected by the finite 
size of the lattices and by the performance of the numer- 
ical algorithms, and it is therefore interesting to compare 
the direct measurement of T with the measurement of m 
in the XY model. In order to measure T in simulations, 
we write it as an integral 



T = / W{z) dz 
10 



(54) 



over an expectation value W(z) 



W{z) 



_ <91ogZ FZS (/3, 7. 
dz 

= 1 P 

Z FZS {p,z) N 3 



E 



<'-M> i ^(^ijkPxjk ~ zm^if ) , (55) 



which can determined with Multicanonical Monte Carlo simulations (for technical details, we refer to Ref. 13). 
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FIG. 5: The tension T at p = 3.0016 (slightly below /3 C , 
i.e. in the Coulomb phase) as a function of N~ 1 . In the 
infinite-volume limit, the tension extrapolates to zero. The 
relationship T cx l/N is equivalent to B tx H 2 predicted in 
Ref.El 



1 .5 



I 

X 


I 1 I 1 I 

: 8 3 


II. 


A 


:16 3 




• 


:24 3 --Jr 




- x / 






I 





, I , I 



3.2 3.4 3.6 3.8 



FIG. 6: The vortex tension T on 8 3 , 16 3 and 24 3 lattices and 
a fit to 24 3 data in the superconducting phase of the FZS. The 
solid line corresponds to the fit in Eq. 102t and the dashed 
line to the XY model mass in Eq. (1441 . 



In Fig. 01 we show the function W(z) measured from a 
16 3 volume as a function of z for several values of (3. 
/3-values in the superconducting phase and close to the 
critical point. 

The vortex tension T vanishes in the Coulomb phase 
and at criticality (3 = f3 c . In Fig. \5\ we show the mea- 
sured value of T as a function of the lattice size N, using 
(3 = 3.0016, which is very close to f3 c in the Coulomb 
phase. The tension decreases with increasing linear size 
N, and assuming finite size corrections of the form N~ a 
we obtain the fit 



T{N) = 0.006(26) + 2.0(7) iV" 



a = 1.1(2), (56) 



with x 2 /d.o.f. = 0.50 for the fit. The data support a van- 
ishing tension in the Coulomb phase close to criticality. 
The value of a is consistent with a = 1 and, at criticality 



and in the dual XY model a = 1 corresponds to a scalar 
correlation length finite size scaling £(N) cx N. Fixing 
a = 1 we obtain the fit 



T{N) 



1.76(3) 
N 



(57) 



at the x 2 /d.o.f. value of 0.65 for the fit. 

In fact, the relation T cx l/N is equivalent to a conjec- 
ture by Son2£ that the flux density B should be propor- 
tional to the square of the external field H. In continuum 
normalization, the flux density corresponding to one vor- 
tex is 

and the external field needed to create it is given by 

where the coupling constant e and the lattice spacing 5x 
are dimensionful quantities. Son defined the constant of 
proportionality C by 

B l^-) 3 CH 2 , (60) 



and argued that it should be a universal quantity. Using 
Eqs. H58[) and (|59|l . we can rephrase this as 



T = C~ 1/2 N- 



(61) 



and Eq. tells us that C = 0.32(1). 

In the superconducting phase, we calculated the vortex 
tension T using 8 3 , 16 3 and 24 3 lattices; the results are 
presented in Fig. as functions of /3. The finite size 
effects between 16 3 and 24 3 lattices are smaller than the 
statistical errors. We fit the 24 3 data in a broad scaling 
interval 3.05 < (3 < 4.1 with the power-law behaviour 

T{[3) = A T ((3-p c yT, (62) 

where vt denotes the tension scaling exponent. Wc 
obtain the parameter values At = 1.24(2) and vr = 
0.672(9) with a x 2 /d-cf. value of 1.13 for the fit. 

The duality (|53|l implies that At and v? ought to be 
equal to the XY model quantities A^ and v in Eq. I|44|l . 
Indeed, we find that v? agrees with the XY model ex- 
ponent v = 0.671, and that At/A^ w 0.94(5). In sum- 
mary, the vortex tension T((3) of the gauge theory agrees 
within two standard deviations of statistical errors with 
the scalar mass m(« = 1//3) of the XY model. This 
constitutes a highly non-trivial test for the methods de- 
veloped in Ref. 20. 



VII. SUSCEPTIBILITIES 

We shall now discuss the susceptibilities Xm 

and xa 

defined in Eqs. l|5rj|) and (|3"7Jl . respectively. We also gen- 
eralize the definition of xa to non-zero momenta by 

/ \ r 33 (p) 

Xa(p) = — , (63) 
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lnx A 




FIG. 7: Logarithm of the gauge field susceptibility lnxA(p) 
[see Eq. 163H 1 as a function of lnp^ 2 in the Coulomb phase 
at P = 2.0 (top curve), at the critical point = j3 c (second 
curve from above), and for two (3- values in the superconduct- 
ing phase. Different symbols correspond to different lattice 
sizes. 

where p\^ 12 ) IS defined as 

2 2 

Pa(i,2) = 2^[l-cosp 4 ] = 2 5^[l-cos(27rJfci/JVi)]. (64) 

i=l i=l 

We use the lowest momentum value for the (1,2) -plane 
momentum, but non-zero values of the z-component of 
the momentum p% = 2irk^,/N^. We expect the momen- 
tum dependence of the susceptibility to be a function of 
the total lattice momentum squared 

3 3 

p 2 A = 2J2[1 ~ cos Pl ] = 2^[1 - coa(2Trki/Ni)}. (65) 

z=l i=l 

As discussed in Sections IIV Al and IIV Bl at zero momen- 
tum xa is finite and non-zero in the superconducting 
phase and diverges in the Coulomb phase, whereas Xm 
vanishes in the superconducting phase and is non-zero in 
the Coulomb phase, being equal to the helicity modulus 
T of the XY theory. 

Our measurements demonstrate the presence of a 
massless photon pole in the Coulomb phase of the gauge 
theory. At f3 = 2.0 - well within the Coulomb phase - we 
determine the susceptibility xa(p) on 16 3 , 32 3 , 48 3 and 
64 3 lattices as a function of small momentum values p. 
The results are shown in Fig. |7|as a function of In p^ 2 . 
The data are fitted with the form 

Xa(p) = ^P- (3 = 2.0, (66) 
Pa 

with Xd o f — 0-67 f° r the fit. This is equivalent to saying 
that in the zero-momentum limit, Xm = 0.46(1) at j3 = 
2.0. This value is slightly smaller than the result from 
free electrodynamics in Eq. I|39|) . 




FIG. 8: The magnetic permeability \ Eq. 

as a 

function of /3 in the Coulomb phase. The curve corresponds 
to the fit in Eq. JB7). 



In fact, since xa(p) is essentially the gauge field cor- 
relation function, albeit in a gauge-invariant form, the 
divergence 1 jp\ demonstrates the presence of a massless 
photon. Thus, we can interpret Xm = T as the corre- 
sponding wave function renormalization Z factor. 

As the critical point is approached from the Coulomb 
phase Xm vanishes with the exponent v, as discussed in 
Section TlV Bl The measured values of Xm are presented 
in Fig. [SI together with the fit to the scaling ansatz 

Xm (/3) = Ar (&-/?)"• (67) 

The parameter values are Ay = 0.46(1) and v = 0.66(2) 
and the fit has a x\ / value of 0.32. By duality, this 
is exactly the scaling of T near the critical point in the 
XY model, and the value of v is indeed entirely consistent 
with arguments that T scales with the exponent vS^ The 
amplitude value A-^ could also be compared directly with 
the XY model, but we are not aware of helicity modulus 
data for the Villain action. 

The anomalous dimension t\a of the gauge field at crit- 
icality can be obtained from the momentum dependence 
of the gauge field susceptibility xa (p) as 

X^pHc^) 1 -^ 72 (3 = f3 c , (68) 

for low p. We show the data together with a power law 
fit in Fig. □ as well as in Fig. ©. The fit to Eq. fBSft 
gives the exponent 

r} A = 0.98(4), (69) 

which is perfectly consistent with the value t\a = 1 pre- 
dicted in Refs. IslIM If we fix t\a = 1, we obtain the 
critical amplitude 

c A = 20.7(3), (70) 
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FIG. 9: The inverse gauge susceptibility \^ (p) at the critical 
point P = 3.002366 for 16 3 , 
lattices as a function of pa 
Eq. fnV 
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The dashed line corresponds to 



FIG. 10: Xa (p) at /3 = 3.02 in the superconducting phase 
for 32 2 x 64, 48 3 and 64 3 lattices as a function of pa . The two 
dotted lines, which are barely distinguishable, correspond to 
fits using Eqs. I|73|l and l|76|l and extrapolate to finite values 
at zero momentum. 



which implies 



xa(p) 



0.0484(7) 



PA 



(71) 



for small momenta. This fit is shown as a straight dashed 
line in Fig. [§] 

In the superconducting phase of the FZS the magnetic 
permeability Xm vanishes at zero momentum, and there- 
fore the photon becomes massive. The gauge field sus- 
ceptibility xa(p) tends to a finite limit as p — > 0. This 
is shown in Fig. where the lowest two data sets of the 
figure correspond to xa in the superconducting phase. 

At the critical point ((3 = [3 C ), xa(p) is proportional to 
pW2-i [ see Eq. (JBSJl]. The phase transition in the FZS 
is continuous, and therefore the analytic form of xa(p) 
at non-zero momenta must interpolate smoothly between 
the critical, massless behaviour and the massive mode in 
the superconducting phase. This reasoning leads to the 
Fisher scaling relation 



1A 



VA 



(72) 



Perhaps the simplest way in which this could happen is 
if the susceptibility has the form 

Xi 1 (p) = CA(^ + n4) 1 - Wa , (73) 

which would correspond to the asymptotic behaviour 

r(r) = r(r,p = 0) - t - va/2 exp(-m 7 r), (74) 

of the photon correlation function. 

In the p = limit Eq. I)73|) relates the photon mass 
and the gauge field susceptibility xa through 



XA 



2-?M ' 
CATTl-y 



(75) 



While the precise form of the ansatz does not affect the 
scaling relation i|72[) . it would lead to a systematic error 
in the determination of m 7 . To estimate these errors, 
we also consider an ansatz based on treating the dual 
XY model as a theory of a free complex scalar field (see 
Appendix j 



n 4 1 
XA{P) = J 

7T CAP 



2p 



— arctan 



(76) 



Again, m 7 gives the exponential decay rate of the cor- 
relator and at the critical point, the ansatz agrees with 
Eq. |[BH||. However, Eq. J7SJ) becomes 

XA=£ — • (77) 

Therefore, we can estimate that the systematic error in 
the determination of m 1 is roughly 3ir/8 — 1 m 20%. 

As an example, we show x~a a t /3 = 3.02 in Fig. 1101 
together with fits of momenta < p\ < 0.5 to Eqs. l (75)l 
and H76|). There is practically no difference between the 
fits, but the fit parameters ca and m 7 are different. 

Performing the fit to Eqs. (|73l) and H7bj) at several (3 > 
(3 C we obtain xa = Xa(p ~^ 0) and m 7 as functions of (3. 
The results for xa(P) and m 7 (/3) are shown in Figs. ITT1 
andHH 

The gauge field susceptibility Xa{P) clearly has a 
power law divergence as we approach the critical point. 
Fitting the ansatz 



X A = A x {(3-(3 c y 



-7A 



(78) 



to the datapoints shown by full symbols in Fig. ^2 we 
obtain 



7A = 0.68(3) for Eq. O, 
j A = 0.70(4) for Eq. tfQ, 



(79) 
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FIG. 11: lnxA as a function of ln(/3 — f3 c ) in the supercon- 
ducting phase, together with power law fits to Eq. 1781 . Here, 
as well as in the other figures, the data points that were in- 
cluded in the fit are indicated by filled symbols and the others 
by open symbols. 
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FIG. 12: The logarithm of the photon (pole) mass m 7 , to- 
gether with power-law fits (dashed lines) and the theoretical 
expectation m 7 = 2T. The circles and triangles correspond to 
zero-momentum extrapolations using ansatze l|73jl and 
respectively. The dashed lines show fits to the filled symbols 
with the XY model exponent v' = vxy- The dotted line is 
the same as the solid line in Fig. 1171 

both of which are compatible with the result ja — vxy 
obtained in Ref. 14. Assuming 7^4 = ^xy, we obtain 

A x = 0.00871(15) for Eq. O, 

A x = 0.00888(17) for Eq. (80) 

The critical behaviour of the photon mass m 7 (/3) (or 
equivalently, the penetration depth) is more subtle. In- 
serting the measured values for j]a and 7^4 in Eq. (|72H . 
we obtain the critical exponent 

v = 0.67(4) for Eq. 



v' = 0.69(5) for Eq. (81) 
This result is fully compatible with the prediction v' — 

8.14.15.26 

Similarly, the scaling ansatze fix the critical scaling 
function m 7 = A 1 (j3 - j3 c y' through Eqs. (O and J77|, 

A 1 = \ - 4.69(8), for Eq. O, 

^ = (l\A = 3 - 94 ( 6 )' for Eq - (82) 

The corresponding critical scaling functions for m 7 are 
shown in Fig. El as dashed lines. The scaling functions 
only agree with direct mass measurements in a narrow 
region near the critical point. In fact, the direct mass 
measurements would favour slightly critical exponents for 
the photon mass 

v' = 0.58(3), for Eq. O, 

v' = 0.62(3), for Eq. (JSJ. (83) 

The reason for this discrepancy is that ca changes rather 
rapidly as a function of (3 when the critical point is ap- 
proached. This is shown in Fig. (|13|) . Using the critical 
value of ca in Eqs. (|75|) and l(77)l is therefore justified 
only very close to the critical point. 

Since the photon mass and the vortex tension scale 
with the same exponent, it makes sense to calculate the 
the amplitude ratio A 1 jAr- We find 

A^/At = 3.8, for Eq. lf?3|). 

A^/At = 3.2, for Eq. fTS). (84) 

These values are significantly greater than the expected 
value of 2, which follows from the assumption that if 
to 7 > 2T, the photon should be able to decay into two 
vortices, 1 and this would result in exponential decay with 
rate 2T. We see two plausible reasons for this behaviour: 
Firstly, it is likely that our ansatze l|73() and (|76|l are 
not of the right form. In this case, it would be im- 
portant to find a theoretically better motivated ansatz, 
which could then be tested numerically by fitting it to 
our measurements. The correct functional form should 
yield Ay/ At = 2. Secondly, the volumes available for 
our analysis may be too small. Vortex-vortex interac- 
tions are presumably not negligible if the lattice size is of 
order l/m 7 , which modifies the "free vortex" behaviour 
m 7 = 2T. 

On the other hand, it may also be possible that the 
system has a massive photon state, similar to resonances 
in particle physics, which would eventually decay into 
two vortices at distance longer than what we can probe 
in our simulations. 

VIII. PHOTON MASS 

In this section, we shall discuss the determination of 
the penetration depth A, or equivalently the photon mass 



13 




10* 



10' 



10 u 



r(r) 




FIG. 13: The value of the constant ca determined from fits 
to Eq. ( 1751 [circles] and Eq. [triangles]. 



FIG. 14: The zero-momentum photon correlation function 
r(r) at j3 = 2.0 in the Coulomb phase. The fitted curves 
correspond to a massless photon. 



m-y = 1/A, directly from the decay of the correlator 
r(r, p) in Eq. (|33|l . This approach is less likely to be 
sensitive to systematic errors than the pole mass deter- 
mination. 

We note that the current operator A*y, i.e., the dual 
of the photon, has odd parity with respect to reflections 
of the gradient angle x +i — 8 X . The more standard com- 
pact functions of the gradient angle 



O (x,«) = if3 'sml 



Oe(x, i) 



-If 

1, 



-0 X ) 

-0 X ) 



(85) 
(86) 



have odd and even parity respectively. The correlation 
functions of these observables map to correlation func- 
tions of 

(x,ij) = e~S sinh(/3P x ,^) (87) 
O e {x,ij) = cosh(/3P x , y ) (88) 

in the FZS, again with definite parity. 

We measure the photon mass from symmetric iV 3 and 
elongated N 3 > N 2 = Ni periodic lattices, using the 
lowest non-zero value of the momentum 



2tt 



(89) 



The asymptotic decay of the photon correlation function 
r(T,p) ~ e - E (p) T i s governed by the lattice dispersion 
relation 



E(P) 



PA(l,2) 



S(P), 



where 



= ^mh(^). 



(90) 



(91) 



The photon self energy S(p) is expected to be small for 
large lattice sizes, and for the data used in this work it 
is unobservable. 




FIG. 15: Photon energy E, defined in Eq. ffity. at j3 = 
j3 c . The straight line is a fit to Eq. 194H and extrapolates to 
vanishing photon mass at infinite volume. 



The photon is massless in the Coulomb phase of the 
gauge theory. In Fig. 1141 we show the correlation func- 
tions r(r) at j3 = 2 measured from 16 3 - 64 3 lattices. 
The curves in the figure are fits to the data, assuming 
vanishing photon mass and no anomalous dimension, 



r(r) = A 



= -r % /2[l-cos(27r/iVi)] 



+ e -(N 3 -T)y/2[l-c OB {2*/N 1 )]'\ + const ( 9 2) 

which yields a perfect fit to the data. 

The photon mass also vanishes at the critical point 
(3 = Pe, but there one must take into account the non- 
zero anomalous dimension t)a = 1. Instead of a pure 
exponential fit, we assume the asymptotic behaviour in 
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FIG. 16: (a) Correlation functions in the FZS and (b) their dual counterpart in the XY model. The data is at /3 = 3.01, in 
the superconducting phase of the FZS. From bottom to top, the curves are the photon, O and O e in the FZS and their duals 
in the XY model. 



Eq. 174(1 and use the ansatz 

r(r) = A[r- l / 2 e- TE + {N 3 -T)- l ' 2 e-^-^ E ] + const 

(93) 

in our fits. 

In Fig. ^] we show the energy of the photon state with 
|p| = 2ir/N, measured from cubical N 3 lattices, as a 
function of l/N. Fitting the data with the form 



E(N) = E(N = oo) + A E N~ 



(94) 



we obtain a vanishing photon energy E(N = oo) = 
0.01(2) and A E = 6.3(6) with x 2 /d.o.f. val ue of 1.20 for 
the fit. The value of A E agrees very well with the value 

27T expected if E(p) = ^Jp\^ 2 y 

Let us now turn to the photon mass in the super- 
conducting phase. At fixed (3 — 3.01 we compare the 
finite-momentum photon correlation functions and cor- 
relation functions of the observables O e and in the 
FZS with their dual counterparts in the XY model. The 
simulations are carried out on 64 x 32 2 lattices with high 
statistics of about 10 6 sweeps. Fig. ED displays three 
|p| = 27r/32 correlation functions for (a) the FZS and 
(b) the XY model. The corresponding correlation func- 
tions of the FZS are not identical to the ones in the XY 
model. Because the correlators have been measured be- 
tween planes rather than points, they are more sensitive 
to boundary conditions, and therefore the discrepancy 
may well persist even in the infinite- volume limit. 

Again we use the ansatz in Eq. I|93() to fit the data. 
The fits are done in the r interval 4 < r < 32. For the 
photon mass, we obtain the values m 7 = 0.08(3) in the 
FZS and m A = 0.06(3) for its dual in the XY model. 
The values are consistent with each other and also with 
twice the XY scalar mass value m(/3 = 3.01) « 0.05(1), 
but the statistical errors are large. 

Mass values from the odd parity O correlation func- 
tions have the values 0.10(4) in the FZS and 0.06(3) in 



the XY model and therefore are degenerate with the pho- 
ton mass. For the even parity O e operator, we obtain the 
mass 0.28(2) in the FZS and 0.25(1) in the XY model. 

In order to obtain the critical behaviour of the photon 
mass m 7 in the superconducting phase, we repeat the 
correlation function analysis at several values of (3 > (3 C , 
using lattices of sizes 48 3 and 64 x 32 2 . The final re- 
sults are shown as circles in Fig. 1171 For comparison, the 
figure also shows m 7 determined from the finite momen- 
tum dispersion of the gauge field susceptibility xa(p) in 
Section TVIII Neglecting the points very close to /3 C and 
assuming a power law singular scaling behavior 



m 7 (/3) =A 1 {f3-f3 c ) 1 



(95) 



a fit in the range 3.02 < [3 < 3.06 to m 7 from the 
correlation functions (squares) yields A~ t = 2.6(5) and 
v' = 0.54(6) with x 2 / d - oi - = 1-3 for the fit. This value 
for the exponent v' sw 1/2 is inconsistent with the XY 
value vxy ~ 2/3. It agrees with the mean field value, as 
was predicted in Ref. |J and observed experimentally in 
Ref.OU 

However, as we already argued in Section IVIII the 
scaling law Eq. 195(1 with v' w 0.5 describes only pre- 
asymptotic scaling, and the true critical exponent v' — 
vxy ■ As the measured values of to 7 are higher than 2T, 
we expect the photon to decay into two vortices, which 
would lead to v' = vxy, but also to to 7 = 2T. This 
behaviour is shown in Fig. 1171 bv the dashed line, and is 
clearly incompatible with the direct mass measurements 
except in very close proximity of the critical point. How- 
ever, the statistical uncertainty of the data in this region 
is too large to justify quantitative comparison. A pos- 
sible source for the overall discrepancy is the fact that 
both the lattice sizes and the distance where the mass is 
extracted from the correlation functions is of order l/m 7 , 
while in order to be able to neglect vortex interactions 
and observe to 7 ~ 2T, distances and lattice sizes much 
larger than this are required. This is very difficult to 
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FIG. 17: The photon mass m 7 in the superconducting phase 
as a function of /3. The squares show the values determined 
from the exponential decay of the photon correlation function, 
and the circles show the pole mass values from fits to Eq. I|7M|I 
in Sect. lVlH The solid line is a power law fit to the measured 
values of m 7 . (Only filled squares are used in the fit). The 
dashed line shows the decay rate m 7 = 2T, which corresponds 
to an unstable photon decaying into two vortices. 

achieve in practice. 

We also remark a further complication in the photon 
mass measurements. Our simulations, which typically 
run for about 10 6 Monte Carlo sweeps, exhibit from time 
to time "exceptional configurations" with large contri- 
butions to the photon correlator at large r-distances. 
We suspect that excitations of the vortex loop network 
are responsible. Indeed, it is easy to see that a vortex- 
antivortex pair wrapping around the finite lattice to the z 
direction also contributes to the photon correlation func- 
tion in Eq. These configurations can significantly 
affect the measurements, unless the lattice size is again 
much larger than the inverse vortex tension. 

IX. CONCLUSION 

In this paper, we have used numerical techniques to ex- 
plore the duality between the three-dimensional integer 
gauge theory and the XY model in the Villain formu- 
lation. Our aim was to identify the ways in which the 
duality manifests itself in the correspondence between 
observables and the critical exponents of the two models. 

First, we used the duality to test the method devel- 
oped in Ref. |2(] for measuring the vortex tension. The 
results agree perfectly with the direct measurements of 
the correlation length in the XY model. Similar methods 
can be used in (and were originally developed for) more 
complicated theories, which do not have an exact dual 
description. 



Second, we investigated the critical behaviour of the 
magnetic field of the gauge theory, which is dual to the 
Noether current of the XY model. We measured the 
anomalous dimension t\a at the critical point, as well 
as the scaling exponents ja and v of the gauge field sus- 
ceptibility and the magnetic permeability. Our results 
are compatible with predictions r\A = ja = and 

These results, together with the Fisher scaling rela- 
tion 17211 , imply that the photon mass (or the penetration 
depth A) scaling exponent v' must be equal to the XY 
model critical exponent: v 1 = vxy& Our results support 
this prediction, but only in extremely close proximity to 
the critical point. 

The non-trivial anomalous dimension makes the de- 
termination of the photon mass extremely difficult and 
prone to systematic errors. The ratio of the photon mass 
to the vortex tension was A 1 jAr ~ 3 — 4 in the vari- 
ous approaches we tried. This means that the photon 
ought to decay into two vortices and the photon corre- 
lator should decay with rate 2T at long distances. Our 
data seems compatible with this very close to the critical 
point, but is not conclusive because of large statistical 
and systematical uncertainties. 

Nonetheless, we believe that the photon does indeed 
decay and the high ratio A 1 jAr is due to an incorrect 
ansatz for the propagator. Were the true form of the 
propagator in the proximity of the critical point known, 
it ought to give A 1 j At = 2, at least if the finite volume 
effects can be avoided. This can be used as a stringent 
test for any theoretical calculation of the propagator near 
the critical point. 

The ratio Ay /At can also be measured in simula- 
tions of the full Ginzburg-Landau theory^ and in prin- 
ciple also in superconductor experiments, since the vor- 
tex tension is related to the critical field strength by 
H c i = T/2ir. We believe that such experiments would 
face the same difficulties in measuring the penetration 
depth as we did with the photon mass, but if a reliable 
measurement can be carried out, the ratio A -1 /H c \ could 
be used to estimate how much closer to the critical point 
one would have to go to see the true scaling in the pen- 
etration depth. As in the experiments in Ref. [ToL the 
penetration depth A = m" 1 we measured in our sim- 
ulations was apparently obeying the mean-field scaling 
v' = 1/2, but it became consistent with the inverted XY 
behaviour A 1 jA T = 2 when {(3 - (3 C )/P C < 0.002. 
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APPENDIX A: CURRENT CORRELATOR IN FREE-FIELD THEORY 

In this appendix we derive Eq. (|76(l by assuming that the vortices are non- interacting free fields at the critical point. 
This means that the dual theory is simply a continuum theory of a free complex scalar field (f> with mass m = T . 
The current operator Jj(x) is defined as ji = lm(j)*di(f>- We want to calculate the current-current correlator 
P)ij(p))- To do this, we write 



(2 ? i+pi)0*(p + q)^(q). 



(Al) 



Then the correlator is 



Ui(p)j j (p')) = lJ ^^(2 ft + P4 )(2r J +p;)(^(p + q)0(q)r(p' + r)0(r)). (A2) 



Assuming that the vortices do not interact, the expectation value factorizes, 

d 3 q d 3 r 



and using the tree-level propagator 



(2tt) 3 (2tt) 3 



(2 9s: + Pi ){2r j + p' j ){4>*{p + q)0(r))(0*(p' + r)^(q)), 



"(pMp')> = (2^) 3 <5(p - P '; 



p 2 + m 2 ' 



we obtain the one-loop integral 



(2q l +p l )(2q j + pj) 



(2tt) 3 (q 2 + m 2 )((p + q) 2 + m 2 ) ' 
Now, this is a tensor that depends only on one vector p, so it must be generally of the form 



(M-P)ii(p)) = A(p)(5 tj B( P ) 



PiPj 



(A3) 



(A4) 



(A5) 



(A6) 



where A(p) and B(p) are functions of the absolute value of p only. The duality in Eq. (|35|l implies that up to a 
constant multiplicative factor and with some constant C, the photon correlator is 



(^(-p^p)) = A(p) + C 



PiPj 
P 2 



B(p) + C 



PrPj 
p2 



Since we know that Bi is sourceless, we must have B(p) — —C, and it is straightforward to calculate its value 

n( \ - EiEl/- ( ynu 1 f d3q (2p ■ q + p 2 ) 2 _ _ m_ 

W " p 2 m P>J J [P>> ~ 4p 2 J (2tt) 3 ((z 2 -|-m 2 )((p + q) 2 +m 2 ) 8tt " 



This means that the gauge field susceptibility xa(p) defined in Eq. 



is 



Xa(p) 



pI+pI 



(B 3 (-p)B 3 (p)) 



Pi + P 2 y 



A(p) 



8tt 



l r 

"^2 



A( P ) 



8tt 



(A7) 



(A8) 



(A9) 



We can calculate A(p) by contracting (ji(-p)jj(p)) and calculating the integral, we find 



with the transverse projection operator, 



(D l)A{p) 



d 3 q 



p 2 q 2 - (p ■ q) 2 



p 2 J (2tt) 3 (q 2 +m 2 )((p + q) 



(10) 



A(p) = - 

and consequently 
1 



m P 4to p 

1 arctan , (11) 

16tt 32tt P 2m' v ; 



xa(p) = 



p — Am p 

— — - m+ arctan— , (12) 

lD7rjr |_ Lp Am 
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up to an overall constant factor, which we parameterize is fixed to m 7 = 2m = 27". 
by ca in Eq. I|7t)[) . The value of m 7 , on the other hand, 
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